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Guided by the analogy to the Bose-Einstein condensation 
of the ideal Bose gas (IBG) we propose a new model for the 
A-transition of liquid helium. Deviating from the IBG our 
model uses phase ordered and localized single-particle func- 
tions. This means that finite groups of particles are assumed 
to be phase-locked. These phase correlations can be related to 
the singularity at the transition point and to the occurrence 
of the superfluid density. 

The model leads to the following results: 

1. A possible explanation of the logarithmic singularity of 
the specific heat. 

2. A characteristic functional form for the superfiuid den- 
sity which yields excellent fits to the experimental data. 

3. A quantitative prediction of a small but nonzero en- 
tropy content of the superfluid component. 



I. INTRODUCTION 



In 1938 LondonS presented his famous suggestion 
about the connection between the A-transition in 
and the Bose-Einstein condensation in the ideal Bose gasci 
(IBG). About ten years later this point of view was con- 
firmed experimentally by showing the absence of such 
a transition in ^He. Theoretically it was supported by 
FeynmanEI who argued that London's view was essentially 
correct. Furthermore, the microscopic IBG provides a 
means of understanding some properties of He II. The 
actual properties of helium are, however, in many re- 
spects quite different from the IBG. In view of this we 
propose a model which is a modification of the IBG in- 
vented such that relevant equilibrium properties of the 
real system are reproduced. 

The IBG can be defined by the microscopic states 



(1) 



which are the symmetrized (S) product states of single- 
particle functions (s.p.f.) (pk of momentum hk. The 
states (0) depend on the coordinates rj (where j ~ 
1, . . . , N) of the atoms and on the parameters Uk — Uk- 
The temperature dependent expectation values (rik) of 
the occupation numbers display the phase transition of 
the IBG. 

Our model is based on two major assumptions: 

1. It starts with an IBG-related guess for the form of 
the many-body states. As in Eq. (P, these states 
depend on occupation numbers rik , but beyond Eq. 
(|l|) they contain correlations due to the use of lo- 
calized and phase ordered single-particle functions. 



2. The phase transition is introduced phenomenolog- 
ically by using the IBG values for (n^). Due to 
this feature our model is closely related to the IBG 
and will therefore be called almost ideal Bose gas 
model, or AIBG. 

The basic idea of the AIBG for introducing the rele- 
vant correlations is the following: Normally exponential 
functions (pk oc exp(iA; • rj) are used for the s.p.f. in Eq. 
(0) . Alternatively one may consider sinus functions which 
depend on the component Xj = x ■ rj of rj as follows: 



^Pk cx sin(9a;j + . 



(2) 



Here i is a unit vector in an arbitrary direction, and 
q = X ■ k. For these s.p.f. we introduce the concept of 
phase ordering (p.o.): 



(phase ordering) . 



(3) 



Phase ordered (p.o.) s.p.f. with the same q (but in general 
different fc) are correlated because their squares \'Pk{'i'j)\'^ 
exhibit common extrema. Thus p.o. leads to an extra 
spatial correlation between all atoms with the same q. 
Such correlations may be important if many atoms share 
the same momentum. Therefore p.o. is a potentially de- 
cisive correlation near the transition point in an IBG-like 
model. 

A p.o. may be introduced by physical (or Dirichlet) 
boundary conditions at the walls of the macroscopic vol- 
ume V. Then the resulting effects are surface effects and 
vanish like V^-^^^ for V oo. In contrast to this, we 
obtain finite effects by assuming s.p.f. ipk{l-P-0.) which 
are not only phase ordered (p.o.) but also localized (1.). 
For this purpose the volume V is thought to be divided 
into V/Vq finite boxes of size Vq. The (^9^(1. p.o) are then 
the s.p.f. which are confined to one of these boxes and 
subject to Dirichlet boundary conditions. For these s.p.f. 
we construct a state of the form dll): 



^Mirj,nk,No) ^SYIYI [</Jfc(l.p.o.)]' 



(4) 



vol k 



The product runs over all momenta and all finite vol- 
umes. The states 'I'm depend on the additional param- 
eter iVo which denotes the number of atoms in Vq- For 
A^'o — > oo these states reduce to IBG states. As a surface 
effect of the finite boxes the l.p.o. introduces correlations 



of the order Nr. 



-1/3 



As it stands, the states Wm cannot be used together 
with a realistic Hamiltonian H. Therefore, we multiply 
them by a suitable Jastrow factor. 
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M, 



F : 



(5) 



We use Jastrow functions f{\ri — rj\) which have been 
determinedu by minimizing the energy Eq = {F\H\F). 

We have now sketched the basic idea of the AIBG, 
namely the modification of the IBG due to a phase or- 
dering of the s.p.f. The resulting model is rather close to 
the IBG; in spite of some formal similarity of the under- 
lying states it is not related to the quasi-particle model 
(note that Y^Uk = N for Eq. (|)). In Sec. || the AIBG 
will be defined in detail. This includes the evaluation of 
the energy 



E{nk) = {F1rM\H\F^M) 



(6) 



Evaluated with {nk)iBG this energy yields a logarithmic 
singularity (Sec. III). 

Since our approach is based on the IBG (above point 
2), the critical exponent /3 of the model condensate den- 
sity has the value 1/2. This differs from the critical expo- 
nent ~ 1/3 of superfluid density ps- Therefore, ps can- 
not be identified with the square of the condensate wave 
function. In the AIBG we argue that due to p.o. also 
noncondensed particles contribute to ps', in any IBG-like 
model {P — 1/2) such a contribution will be required in 
order to account for the experimental value i> ~ 1/3. This 
contribution leads to the peculiar consequence that the 
superfluid component has a nonvanishing entropy Sg 0. 
A prediction for Ss will either falsify the model right 
away, or provide a crucial test of the model. Therefore 
we extend the model in Sec. ^ to VI such that ps and Sg 
can be calculated. 



Section [V shows how the localization of the s.p.f. can 
be reconciled with a macroscopic p.o. and introduces 
generalized s.p.f. which are able to describe the coher- 
ent motion of noncondensed particles with the conden- 
sate. Section ^ presents the model expression for ps and 
shows that it yields excellent fits to the experimental 
data. With the model parameters fixed by this fit we 
obtain in Sec. VI the crucial model prediction for Ss ^ 0. 
This prediction is compared to the available experimen- 
tal data; it is found that the prediction is at the border 
of present-day experimental detectability. 



II. THE ALMOST IDEAL BOSE GAS MODEL 
(AIBG) 

For a complet e de finition of our model we defin e the 
l.p.o. s.p.f. (Sec. II A ), calculate the energy E (Sec. OB) 
and specify the assumptions about the expectation values 
of the parameters (Sec. IIC ). 



A. Localized, phase ordered single-particle functions 

The l.p.o. s.p.f. in Eq. (^ are formally defined as fol- 
lows: The macroscopic volume V is divided into V/Vq 



cubic boxes of size Vq. Within each box we define or- 
thonormalized s.p.f. which obey Dirichlet boundary con- 
ditions at the walls: 



V5fc(l.p.0.) = < 



]^sin(fci 



i=l 







for r e Vq , 



elsewhere . 



(7) 



The Cartesian components ki of k are restricted to the 
discrete values g„, 



Qn — n ■ Afc, n = 1, 2, 



AA: = 7^/yo^/^ (8) 



The many-body states of the AIBG are now defined by 
Eq. d) with Eqs. (|) and (0). 

This construction is chosen for simplicity. It contains 
artificial aspects like the cubic shape and the identical 
size of the boxes. The underlying physical picture is that 
there are finite regions (of various shapes and sizes) in 
which the atoms lower their free energy (see Sec. IVA) 
by assuming p.o. in just one direction. For macroscopic 
physical quantities we have to average (explicitly or im- 
plicitly) over such regions, or over the boxes used for Eq. 
(0) . Formally we perform this averaging by replacing the 
discrete momentum sums by integrals: 



E 



(9) 



The justificati on of the lower integral bound will be dis- 
cussed in Sec. IV B. The averaging over boxes of finite 
size and shape restores also the translational and rota- 
tional invariance of the system as a whole. The occupa- 
tion number for ipk will therefore depend only on k but 
not on the direction of k. 



B. Energy 



We evaluate the energy Eq. (0) for the Hamiltonian 



(10) 



with a realistic atom-atom potential u(r), for example of 
Lennard-Jones type. p 

The ground state (g.s.) may be approximatedH by tf'o = 
F. The corresponding g.s. energy Eq can be written as 

EoiV,N) = {F\H\F) = — J d'ru{r)go{r), (11) 

where go{r) is the pair correlation function for and 
{t(r) is given by 



u{r) = u{r) - — ln/(r) , 
2m 



(12) 
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Sensible results for Eq and go {r) are obtainedB by a sim- 
ple ansatz for /(r). 

The following evaluation of the energy E ior >F = F'Pm 
is based on the fact that the main correlations are deter- 
mined by F. The additional correlations due to •Z'm 7^ 1 

— 1/3 

are small, they are of the order A'q <C 1 for the rel- 
evant p.o. effects. The pair correlation function g{r) for 
•Z' = F'Fm can therefore be written as 



gir,nk) ^ goir) + Ag{r,nk), Ag < 1 , 



(13) 



We relate Ag to the pair correlation function gu of tZ/j^. 
(An explicit expression for and details of the following 
derivation are given in Ref. |^). Writing — 1 + A^m 
it follows that Ag for A^m ^ 0, or Ag cx A^m- To 
lowest order in /(r) Eq. (H) yields 



Ag{r, Uk) ~ f{r)^ Aguir, Uk) . 



(14) 



An alternative approximation is Ag ~ go AgM] the dif- 
ference to Eq. (|lj) marks the uncertainty in this approx- 
imation. For our purpose both approximations will turn 
out to be nearly equivalent. The relation between Ag 
and AgM cannot be evaluated exactly; for a discussion of 
this problem see Ref. |^. 

The potential energy for F^m is Epot = (N'^/V) ■ 
Jd^rgu. The kinetic energy can be divided into two 
parts: The Jastrow factors are taken into account by 
replacing u in -Epot by w, as in Eq. (pT|). The kinetic 
contributions due to •Z'm 7^ 1 lead to ^ Sk nk where 
Ek = h^k'^/2m. Here m is an effective mass (Ref. |^) 
but for simplicity we do not introduce a new symbol. We 
emphasize that m 7^ for fc 0. This can be related 
to the fact that the static structure function S{k,T) is 
finite for fc — > and T ~T\. The complete result reads 



E{V,N,nk) = {F^m\H\F^m) 

= ^ j d^r g{r, Uk) u{r) -f ^ Efc rife 
= Eo{V,N)+EM{V,N,nk). (15) 



The energy Eyi — E ~ Eq is obtained by inserting Eqs. 
([|), (0) and dnl): 



EM{V,N,nk) ~ ^Efe nfc 



(16) 



+ ^ j d\ f{rf AgM{r,rn,)u{r). 

The decisive contributions in A^m are long-ranged be- 
cause the low momenta dominate near the transition. For 
the intended evaluation of the critical behaviour we may 
therefore use the approximation 



Agyiir., Uk) ~ A5m(0, Uk) 



(17) 



EM{V,N,nk)=Y.eknk + ^ — AgMiS),nk), (18) 
^ — ' 2 V 

k 

where v = V/N and 



2E, 



g.s. 



= -14.3 /crK. 



(19) 



Using McMillan's parametersB, the first integral yields 
— 14.8A:bK- Within the uncertainty of Eq. ( p^ it may 
be replaced by the second integral. This integral is the 
theoretical g.s. energy ([ll|). By equating this energy to 
the experimental ground state energy £'g g wc fix the 
value of wq/v. 

Up to a constant, the energy Em (^ can be expressed 

by 



E, 



M 



{>Fm\H\'Fm) with 



'^-^ 2m 



(20) 



This means that within sensible approximations the en- 
ergy Em in E — Eq + Em can be calculated as the expec- 
tation value of a simple product state S'm with a simple 
Hamiltonian Hm- It should be noted that Hm is solely 
introduced for evaluating Em , it cannot serve as a model 
Hamiltonian in other respects. In particular the stability 
of the system is guaranteed by the part Eq of the to- 
tal energy E = Eq + Em ', the corresponding part of the 
Hamiltonian is not contained in Hm- 

The attractive zero-range force in Hm has the following 
meaning: The l.p.o. implies an enhanced probability of 
finding two atoms near together. Due to the Jastrow 
factors in Eq. (^) this enhancement is effective only in 
the range where the realistic interaction is attractive; the 
strength of this attraction is measured hy wq/v. A zero- 
range force can be used because for the enhancement 
is approximately r-independent up to a few A. 

The evaluation of Eq. ( pO| ) is straightforward. The two- 
body matrix elements Wkq of the (5-force with the s.p.f. 
(0) are given by 



Wkq 



{ipkiPq\wQ6{ri - r2)\ipk'Pq) 



Wq 

Vo 

Wq 

Vo 



n 

1=1 



1 + 2 ^i''<]' 



1 ^ 



(21) 



This yields 



The somewhat artificial construction Eq. (0) implies p.o. 
in 3 orthogonal directions. In contrast to this we as- 
sume that the p.o. is actually realized in just one direc- 
tion. Therefore we omit products of (5-functions in the 
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last expression of Eq. (^Tj) which originate from simulta- 
neous p.o. in 2 or 3 directions. The remaining (5-terms 
describe the mutual correlation of all particle with the 
same momentum q in one direction. The number Vq of 
these particles is 



= ^ rtfc where k = \k\ = ^ + + k^ . (22) 



The obvious form of the correlation energy is -Econ- oc 
] v^. The actual result follows from Eq. ( po|) with Eqs. 



(^) and (12^): 



3 Wo 1 



2v M 



1/3 




2 no 1^0 



(23) 



For the replacement (^) of the sum by an integral we took 
into account a possible finite condensate fraction uq/N. 
The index zero of uq and vq stands for the lowest possible 
momentum value. 

The complete expression for the energy E of Eq. ( |l5| ) 
becomes 

EiV, N, nu) =Eo+Em=Eo+ Eibg + ^corr ■ (24) 
Here Eq — Eo{V,N) is given by Eq. (|lll ) and Eibg by 

(25) 



EiBG{V,N,nk) = £fc nfc . 



The critical contribution in Eq. ( [2^ ) is the correlation 
energy Ecorr (23). In the evaluation of Eq. (20) leading 
to Eq. (|2^) we dismissed a constant term stemming from 

the constant in Eq. (pT|), the higher order terms (oc Nq^^^ 
or Nq^) because they come from p.o. in more than one 
direction, and a term oc rig because it is noncritical (that 
means eventually negligible near the transition point). 
For the finite boxes the momentum sums have to be 

1 /3 

evaluated with Afc = tt/Vq ; the subsequent summa- 
tion over all boxes yields, however, a factor V/Vq. For 
plain momentum sums (without Jfe^g^ -terms) this proce- 
dure can be abbreviated by performing the common sum- 
mation with Afc = tt/V^/^. The final momentum sums 
ins Eqs. ( p3| ) and (|2^) are performed in this common 
way. This implies that the occupation numbers nk do 
no longer refer to a finite box (as in Eq. (^) but to the 
macroscopic system (with V and N). In Eq. ( p3| ) this 
implies ^ oc N^l^ and an A^-independent Ecoi-r/N. 



C. Statistical assumptions 

In order to calculate the thermodynamic energy 
E(T, V, N) from Eq. ( |2^ ) we need the temperature de- 
pendent expectation values of the parameters of W. Our 
model uses the IBG expression for (uk) and treats A^'o 



as an adjustable constant. We discuss and specify these 
assumptions. 

Any model of the A-transition except the IBG intro- 
duces the phase transition phcnomenologically. In the 
AIBG this phenomenological assumption is the use of 
the expectation values (rifc)iBG- This assumption can be 
made plausible (to some extent) by observing that the 
additional contribution E 
small in the sense 



Er. 



Eibg 



0{y) where y 



\w,lv\N-''^ 



< 1 . (26) 



The actual parameter values will yield y = 0.13, Eq. (|42 
For A'o ^ oo or y ^ one obtains Em Eibg and thus 
{nk)iBG- (The nfc-independent contribution Eo(V, A^) in 
E is without influence on the expectation values) . There- 
fore, we expect (uk) = {nk)iBG [1 + C'(?/)]. The evalua- 
tion of Ecorr with (nfc)iBG IS then valid up to first order 
in y; it is this contribution which yields the logarithmic 
singularity. Another argument for using the (nfc)iBG is 
that the additional term ^ correlates many s.p. states 
with many others, and that it will be without much in- 
fluence on the occupation of a particular s.p. state. 

The AIBG assumes expectation values (uk) of the IBG- 
form, 



1 



[rik) 



exp [{sk - fJ-)/kBT] - 1 
1 

exp (x^ + T^) — 1 



(27) 



Here ^ is the chemical potential and fca is Boltzmann's 
constant. We have introduced the dimcnsionlcss quanti- 
ties t"^ = —^/k^T and 



A|fc| 



with A = 



27rft 



v27rTOfc^T 



(28) 



The transition temperature of the IBG is given by the 
following condition for the thermal wave length A — A(T): 



A(r,) = [<(3/2)]i/3 



(29) 



where C(3/2) = 2.6124 denotes Riemann's zeta function. 
In applying the AIBG to the real system we will identify 
Tx with the actual transition temperature; formally this 
can be achieved by inserting a suitable effective mass 
TO — m{v) in A of Eq. (p8|), Ref. ^. In the following we 
use the relative temperature 



t 



T-Tx 



(30) 



For t > the chemical potential is determined by the 
particle number condition = A^. For t ^ this 

condition yields ii — Eq —> and thus no — *■ oo. For 
simplicity we rename /i — £o by /i; it is this new /i which 
vanishes at t = also for Eq 0. 
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We are interested in the critical properties. Therefore 
we expand fj, or, equivalently, r for |t| <ti 1: 



at + ht^ + . 

a'\t\ + b'e + 



(i < 0) . 



(31) 



This expansion is directly related to the model conden- 
sate fraction 



N 



= 1 - 



nk) 







f\t\+gt' + 



{t > 0), 
{t < 0). 



(32) 



The coefficients /, g, . ■ . are determined by a', in 
particular / = 3/2+20Fa7C(3/2). The IBG yields finite 
values for a, b, . . . , and a' — b' = . . . = 0. 

Deviating from the IBG, the AIBG admits a coefficient 
a' 7^ in Eq. (|l|) or / > 3/2 in Eq. (H). This devia- 
tion is introduced phenomenologically for the following 
reason: For t > the evaluation of ^ with (nfc)iBG 
yields straightforwardly a logarithmic singularity. This 
singularity can be directly continued to t < by admit- 
ting a' 0; the singularity itself does not depend on the 
value for a' . The form of the expansion ( ^ ) is to some 
extent plausible because of its simplicity and symmetry. 
Physically a' 7^ means an energy gap of a'^t^k^T for 
the condensed s.p. state for i < 0. 

For evaluating the thermodynamic energy we need also 
the expectation values of products of rifc's. (We do not 
need, however, the full statistical information contained 
in the density matrix). For such products the IBG yields 



(nfejiq) = (nfc)(nq) {k ^ q) 



(33) 



Compared to the IBG, the AIBG states depend on one 
additional parameter, namely N^. In principle, this pa- 
rameter (as well as the occupation numbers) should be 
determined from the condition of minimal free energy 
F yielding a temperature dependent expectation value 
(iVo). Instead of this, we treat iV p as a n adjustable con- 
stant. Simplified estimates (Sec. [VA and Ref. ^ yield 
a sensible value for A'o at i = 0. 



III. LOGARITHMIC SINGULARITY 

We evaluate the thermodynamic energy E{T, V, N) = 
{E{V, N,nk)). According to Eq. (|^ it is of the form 

EiT,V,N) = EoiV,N) + EiBciTXN) + E,,„{T,V,N). 

(34) 

The term Eq (V, N) is essential for a realistic compress- 
ibility but it does not contribute to the specific heat. The 
energy Eyqg = J2^k{nk) is calculated as in the IBG. 
Equations (|2^) and ( p3| ) yield for the decisive contribu- 
tion Econ, 



N 



2v ^l/3 7V5/3 



2(no)(t^o) 



(35) 



Leaving away unessential constants we sketch how i?corr 
yields a logarithmic singularity: For |t| — s- and x — s- 
the occupation numbers ( p7| ) behave like (rix) ^ l/[x'^ + 
t^] and (z/,) - Jdkk{P + q^+e)-^ ~ ln(g2 +i2). It fol- 
lows that {no){vo) ^ \t\\n\t\ and ^(j^,)^ /dg[ln(g2 + 
t^)Y ^ const. + |i|ln|t|. 

We present now the results of the detailed calculation. 
The {vq) of Eq. ( ^ with Eq. (^) can be evaluated an- 
alytically. 



In [1 



exp I 



-X 



(36) 
and 



We used the dimensionless quantities r of Eq. (|3]_ 
X = gA/Vivr. The integral in Eq. ( p5| ) becomes 

^|:W^^-|CW2)-»'3(0'7W, (37) 
where i(0) = 8.30 and 

L(t) = y da;(^ln[l-exp(-a;2 -r^)] j (38) 
= L(0) + 47rTlnr -|- 47rT [ln(2) - l] + ©(t^) . 



The expectation value (z/q) is given by Eq. (|3^) with x = 
because eo has been absorbed in = —{fi — £Q)/k^T 
(remark after Eq. (|30|)). Using {vq) ^ In |t| (np) ^ \t\ 
and T ^ \t\ we see that both terms in Eq. (pq) yield a 
contribution ~ t • In The critical {\t\ ^ 1) behaviour 
of the specific heat reads therefore 



(39) 



_ J -A-\n\t\+B + .. 
~\-A'- In |t| + B' + . 

The detailed evaluation of Eq. (|3^) yields 



(t > 0) , 
(i < 0) , 



A 



9C(3/2)-2/3^ 



(40) 



The terms depending on the parameter a' cancel. We 
insert \wo\/v of Eq. (^ and Tx = 2.17 K in the result 
(!)• 

Experimentally a logarithmic singularity has. been 
measured for cp over several decades by Ahleraj. For 
saturated vapour pressure we may use cy ~ cp. (The 
consistency of a logarithmic form for both, cy and cp, 
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has been demonstrated by Lee and Pufffl). The com- 
parison of our result (^ ) with the experimental values 
A A' 0.63 fcs (Eq. (48) of Ref. |) determines the 
model parameter iVo, 



1/3 



50 (adjusted) 



(41) 



This fixes also the small parameter y of Eq. which 
is a quantitative measure for the deviations of the AIBG 
from the IBG, 



'1/3 



0.13. 



(42) 



Experimentally_one finds a small difference between A 
and A' (Ahlerstl reports 5%) and an indication of a cut 
of the logarithmic singularity (slightly negative critical 
exponentscl for cp). In the AIBG a difference between A 
and A' could be obtained if corrections to approximation 
( p7| ) are taken into account which differ for t > and 
t < 0. A possible cut of the singularity is connected with 
the lowest possible values of the momenta. 

The coefficients B and B' in Eq. (39) can be calculated 
from Eq. (|35|), too. They depend on the value of t he pa - 
rameter a' which will be determined only later (Sec. |V Q ). 
Moreov er a generalization of the s.p.f. to be introduced 
in Sec. [V C will yield another contribution to the jump 
in B' - B. 



IV. PHASE ORDERING 

We review various aspects of the p.o. assumed in the 
AIBG. First we sh ow that this p.o. is favoured by the free 
energy (Sec. [V A ). Then we review the necessity of the 
localization of the s.p.f. for o btaining finite correlation 
effects due to p.o. (Sec. IV B| ); this leads to the notion 
that the coherence range of p.o. appro aches infinity at 
the transition point. Finally (Sec. IV C| ) we connect the 
p.o. with a potential superfluid flow. 



N, 



A. Free energy due to phase ordering 



The l.p.o. leads to a lower energy, Econ/N ^ {wo/v)- 

1 /3 

< 0. At the same time such an ordering reduces 



the entropy, AS'(p.o.) < 0. We show that AS'(p.o.) is 
relatively small and that p.o. is therefore favoured by the 
free energy. 

Arbitrary phases in Eq. (|^) are statistically equivalent 
to an equal weight of the two independent choices for 
the phase, 0o as in Eq. (^) or 0o -I- 7r/2. We consider 
p.o. in just one direction. For each k we have then two 
s.p. states. The number of possible distributions of Uk 
bosons on 2 states is Uk + 1. The p.o. requires that 
all Tifc bosons go into the same s.p. state; compared to 
a statistical distribution this implies the entropy change 
As(p.o.) — — /cBlii("-fc + !)• On the other hand, the Uk 



atoms gain by p.o. the energy Ae(p.o.) {Ecoii/N)nk. 
The free energy change A/ is then 



A/(p.o.) ^ - 



vN, 



-1/3 



Uk + ksT \n{nk + 1) . (43) 



In spite of the smallness (^2|) of the considered corre- 
lations the p.o. is favoured (A/(p.o.) < 0) as soon as 
nfc 3> 1. The boson property of the atoms is decisive for 
this conclusion: It implies that |As| is not proportional to 
Uk but only to In ; this behaviour may be paraphrased 
by 'bosons like to go into the same state'. For the rel- 
evant lowest s.p. states {k ~ 0) the condition nfc ^ 1 
is fulfilled for i ^ 1 . Thus p.o. will indeed be adopted 
when the transition point is approached {t 0"*"). 

As a simplification our ansatz (^) assumes p.o. for all 
s.p.f. Summing over all momenta and including numeri- 
cal factors one finds then Econ — TAS{p.o.) for T 
and for the A'o of Eq. (p|). Reversely, the condition 
-Bcorr ~ T\ AS'(p.o.) constitutes a rough theoretical es- 
timate for the model parameter TVqj leading to an TVq of 
the order (Ell). 



B. Phase coherence volume 

In the Introduction we argued that for finite correlation 
effects due to p.o. the s.p.f. have to be localized. We 
review this point and find that the localization is not 
required for the lowest s.p. state. We present then an 
estimate for the extension of this lowest s.p.f.; within 
this extension the p.o. will be coherent. 

— 1/3 

The finiteness of the correlation energy (cx iVg ) is 

1 /3 

due to the finite spacing Ak = t^/Vq of the values in 
the sum i?corr ^ Tli^q ■ Noulocalizcd s.p.f. would lead to 
the same -Econ- provided that only s.p. states with 



Ak, 



(n = 0,l,2,...). 



(44) 



are occupied. In the macroscopic system the possible q- 
values are, however, dense, and the entropy drives the 
particles into the occupation of all available states. This 
is the reason why a finite Afc can be realized only for 
s.p.f. localized within a volume Vq = (ir/Ak)^. 

This argument shows that we may either start from 
finite volumes or from finite Afc, the second feature fol- 
lows from the first one. There is, however, one difference: 



Starting from the finite spacing Afc in Eq. (44), the nat- 
ural condition for the lowest (/-value is qo < Ak. The 
corresponding s.p.f. can then not be localized within Vq; 
as just stated there is also no such need for a localization 
for the required correlation effect. 

In view of this we introduce the following modifica- 
tion of the states (^: Only the s.p.f. with qn>o are lo- 
calized within Vq, the s.p.f. with q^ have a larger ex- 
tension. We present a crude estimate for the volume 
Vc of the lowest s.p.f. ifo which becomes the conden- 
sate state for i < 0: Let Vc be some multiple of Vq, 
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that means Vc = WVo- In a volume Vc there are W 
s.p. states below Afc out of which only one is occupied. 
A redistribution of the no atoms over these W states 
would increase the entropy by As ~ fcBlii(no)'^. At 
the same time these atoms would loose their correla- 

— 1/3 

tion energy, Ae ^ no\wo/v\NQ . The stability con- 
dition Tx As < Ae yields an upper bound for W. Using 
{no{t > 0)) - t-^ and {no{t > 0)) = 0{N) we obtain 



r 



(i > 0) , 
it<o). 



(45) 



This result means the following: Without changing our 
previous calculations we can modify the states (^ in such 
a way that the lowest s.p.f. ipo has the large and even- 
tually infinite volume Vc, whereas the srpl. with > 
are localized. This modification of Eq. (|) has important 
consequences: 

1. Within the volume Vc the direction of p.o. is defined 
by (fio- Therefore, Vc{t) of Eq. (^) is the phase 
coherence volume of the considered p.o. 

2. The p.o. constitutes a symmetry breaking because 
in an infinite system the atoms are free to adopt 
arbitrary phases (pj in Eq. (H). (The average over 
the s ize and shape of finite boxes (assumed in Sec. 
[I A ) does not restore this symmetry). The result 
([4 51) means that this symmetry breaking changes 
its character from local to global at the transi- 
tion point. Approaching the transition point {t — > 
0+) the directions of p.o. of neighbouring Vq's get 
aligned, and for t < the coherence is potentially 
infinite. 

3. Strictly finite and separated boxes imply a momen- 
tum cut at Afc and consequently a cut of the log- 
arithmic singularity. This cut would be described 
by a lower bound Afc in the integral in Eq. (^. For 
Eq. (^4|) we revise the replacement Eq. (|^) of the 
momentum sums by integrals: 




1 
Afc 



dq 



(46) 



The brackets indicate an average over possible go's 
with go l£ ^k. The lowest value for go is deter- 
mined by Eq. (^5|). A lower bound go.min ~ (for 
i > 0) of the integral does not lead to a cut of the 
logarithmic singularity. 



C. Complex phase ordering 

The phenomenon of superfiuidity can be explained by 
assuming a macroscopic wave function ^. The phase ^(r) 
of the complex ip yields the velocity Vs = {h/m)'V<P of 
a potential superfluid flow. In the AIBG the condensed 



particles for a macroscopic wave function for t < 0, too, 
Eq. (^). In order to join the common description we 
replace the real s.p.f. ipo by the complex s.p.f.: 



(fio = (fio exp [i^(r)] . 



(47) 



This replacement does not affect the previous calcula- 
tions; moreover, the Jastrow factors in Eq. ^ arc with- 
out infiuence on a flow due to V<P ^ 0. With Eq. (|4^) 
all results obtained from the assumption of a macroscopic 
wave function apply to our model, too. Besides this basic 
conformity the AIBG leads, however, also to a peculiar 
modification which will be discussed in the following. 

In an IBG-like model (with (uq) ~ \t\ or /3 = 1/2) we 
have to assume that noncondensed particles contribute 
to the superfluid density ps in order to get agreement 
with the experiment {ps ^ jip/'^ or v = 1/3). The con- 
sidered p.o. should in some way imply that noncondensed 
particles move coherently with the condensate. We intro- 
duce the possibility of a net current of the noncondensed 
particles by replacing the original real s.p.f. ipk, by 



Vk = exp [i<Pk{r)] . 



(48) 



A coherent motion with the condensate can now be de- 
scribed by the condition = <1>. We call this condition 
complex phase ordering (c.p.oj named after the com- 
plex phase factors in Eqs. ( ^Tj) and (|4^). In contrast 
to this the p.o. (^ will be called real phase ordering 
(r.p.o.). We note two points: The spatial correlations 
of the s.p.f. are unchanged by the replacement (^7|), ( |4^ ) 
because = Secondly, for a coherent macroscopic 
flow it is sufficient that the field <P{r) is macroscopic; it 
is then not necessary that all contributing s.p.f. are of 
macroscopic range. 

The replacement (^^ is sug gest ed by the familiar pic- 
ture of a superfluid, and Eq. (|4§| ) with <Pk = <P is nec- 
essary to reconcile (3—1/2 (model) with = 1/3 (ex- 
periment). The fairly straightforward introduction of the 
phase factors (|4^) and ( p7| ) determines the crucial model 
prediction for Ss ^ 0. The fit to the experimental Ps 
determines the extent of the c.p.o. and consequently Sg- 

Presenting a peculiar prediction like Sg we feel 
obliged to make it quantitative. For this purpose we have 
to write down explicit expressions for ps and Sg. This 
requires additional assumptions which are introduced by 
plausibility arguments. As explained in the last section, 
the quantitative result for Ss will not depend sensitively 
on these assumptions because the model expression for 
Ps is fitted to the experiment. 

We require that the new modes 'Pk do not destroy the 
correlation energy due to the r.p.o.; this determines the 
extent of the c.p.o. The field <P{r) contains degrees of 
freedom which will be thermally excited. Let kc{t) be 
the average amount of the momenta of this field 



kc{t) ^ {\V<P{r)\) . 



(49) 



Because of = \(p\'^ the fields <P and <Pk do not directly 
infiuence the spatial correlations. However, for = 
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the noncondensed states k with < h'^k^/2m would 
have a lower energy than the condensate state. This 
would invalidate the assumption about the occupation 
pattern and destroy the correlation energy. On the other 
hand, for k ^ kc the ^fc will have only minor influence 
on the occupation pattern. Simplifying and quantifying 
this qualitative argument leads to the following condition 
for c.p.o.: 



^r) {k < k,) , 
arbitrary (fc > kc) . 



(c.p.o.) (50) 



This condition means that the localized s.p.f. with k < k^ 
adopt within their range the (macroscopically defined) 
phase 'P of the condensate. The physical reason for this 
c.p.o. is that the correlations due to the r.p.o. (favoured 
by the free energy) must not be destroyed; in this way 
the real and complex p.o. are connected to each other. 



Using Eq. ( p3| ) and ( V^th ) = we evaluate the statis- 
tical expectation value of Eq. (^). Writing ( jg) = p^Vs 
with Vg — (?i/m)V<?s we obtain for ps'. 



7 = ^'^"°^ 



Po 



(54) 



This is the AIBG expression for the superfluid density. 
In a number of points we summarize its implications: 

1. The density ps is composed by all atoms adopt- 
ing the phase <P{r). It consists of the condensate 
density po and the comoving density pc- The tem- 
perature dependence of this composition (Fig. ^ is 
obtained by fitting Eq. (jsj) to the data. The co- 
moving density pc has an internal structure which 
implies a nonvanishing entropy Sg ^ of the super- 
fluid component. 



V. SUPERFLUID DENSITY 



A. AIBG expression 

Using the new s.p.f. (^^ and ( ^ ) the states (jHj) become 
for t < 0: 

tf' = F5((po exp[i<?(r)])"" [] (^^ exp [i<?fe(r)])"\ 



2.0 1.9 1.8 T{K) 1.7 



fc,vol 



(51) 



For these states we evaluate the quantum mechanical ex- 
pectation value of the current operator. Taking into ac- 
count the c.p.o. (|50|) we obtain 



jir,nk) = ( !f 



N 



-^V,,5(r-r,) 



+ c.c. 



j=i 



(52) 



p h 
Nm 



no 



All derivatives of the real functions {F, (po and ipk) 
are cancelled by the complex conjugate (c.c.) term; 
the only surviving contributions are the phase deriva- 
tives. Each such derivative is accompanied by a factor 
{'P\S{r — rj)\<P) = p{r)/{Nm). The mass density p is 
assumed to be a constant in the following. 

The contributions of the nonordered phases are aver- 
aged out statistically. This is also the case for the ther- 
mal excitations of the <?-field, (V<?th) = 0. A coherent 
nonvanishing flow can only be obtained by an additional 
nonequilibrium contribution <Ps to ^, 



<P = <P 



'th 



(53) 



The flux due— to <?s might be relatively stable if it is 
small enoughlij. This implies |V<?s| <C fee where fee = 
(|V<?th|). 




FIG. 1. Composition of the superfluid density ps/p (solid 
line) according to Eq. (^^. The ratio ps/p is the sum of 
the model condensate po/p (short dashes) and the coherently 
comoving density pc/p (dash-dotted line). Also shown is the 
simplest 1-parameter fit ps/p = ai \t\'^^^ (long dashes). In the 
given scale the fitted pa/p (solid line) coincides with the data. 

2. Expression ( p4| ) is the simplest possibility to ac- 
count in an IBG-likc model with (rto) ~ \t\ for 
Ps ^ Fitting Eq. ( p^ ) to the experimental ps 
fixes fcc(i), and consequently the difference between 
the densities po and Ps and the model prediction for 

3. The reproduction of the data with Eq. ( p4| ) requires 
asymptotically kc{t) ^ \t\'^^^. Theoretic ally t his 
asymptotic form is made plausible in Sec. 



IVB. It 



leads to a certain functional form of ps which yields 
excellent fits to the experimental data (Sec. 



VC) 



4. The result (^ contains a specific picture for the 
relation between the condensate fraction ric and the 
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superfluid fraction Ps/p- The explanation of super- 
fluidity is based on the assumption of a macroscopic 
wave function. This connection implies, however, 
an open question: What is the quantitative rela- 
tion between nc(T) and ps{T)/p, and in particular 
the relation between the T — values rtc ~ 0.1 
(see Ref. |l|) and ps/p =1? The AIBG pro- 
poses the following picture: The model condensate 
n-c,M = {no)/N — po/p is the fundamental quan- 
tity. It is depleted by the Jastrow factor F in Eq. 
( pT[ ) to the real condensate ric, for T = the value 
n-cM = 1 is reduced to ric ~ 0.1. The factor F 
in Eq. ( |5l] ) does, however, not deplete the current 
due to phase factors. Therefore n^M 1 implies 
Ps/p 1. On the other hand, for \t\ <C 1 the 
contribution pc ^ po dominates ps- 



B. Effective Ginzburg-Landau model 

Ginzburg and SobyaninEl have proposed an effective 
Ginzburg-Landau functional for the free energy _Fgl as 
a function of the order parameter field ip. In this ansatz 
singular coefficients (like 1^1^^"^ for the j^/ip term) are in- 
troduced in order to yield the right critical exponents. 
We follow this kind of approach for investigating the re- 
lation between the critical exponents of the order param- 
eter and the superfluid density. This detour is, however, 
not required for the final result for Sg (see point 2 of Sec. 
VA| ). 

For discussing the fluctuations of the <^-field we define 
the order parameter field by 



exp [if?(r)] 



(55) 



For a definite phase the state .-(eTI) must be replaced by 
the appropriate coherent stateli37^This means that uq in 
Eq. ( p5| ) has to be understood as a quantum mechanical 
expectation value of the occupation number no in such a 
coherent state. 

The statistical expectation value (ng) ~ |<| can be 
obtained by minimizing the common Landau energy 
Fl/V — Rt \ip\'^ + U \ip\'^ (with regular coeflicients R and 
U). Adding a naive kinetic energy term (?i^/2m)| VV'P 
leads to a violation of scaling invariance and to wrong 
critical exponents. The decisive feature of the AIBG ex- 
pression (^4|) is that the mass density comoving with 
is Po + Pc instead of po = to( IV'P) alone. This suggests 
the following effective Ginzburg-Landau ansatz 

E5^^^PL\Vi,?+Rt\^' + U\M\ (56) 
V 2m Po ^ 

Leaving away unessential constants, the asymptotic be- 
haviour of the mass density pc in Eq. (54) is 



Jk<ka 



(57) 



This, p ~ Pc for 1*1 and Ec _ 
asymptotic kinetic energy in Eq. (pq), 

^ iWiA') = ^ {\v<P\') - k! . 

Po ' m ^ ' 



(Eoh determine the 



(58) 



Scahng this kinetic part of Fql with (Fl) ~ yields 



1^1 



2/3 



(59) 



This implies ps ~ Pc ~ \t\^^^ and a singular mass coeffi- 
cient ps/po ~ \t\~'^^^ in Eq. (H). Therefore, Eq. @) is 
an effective Ginzburg-Landau ansatz in the sarne-sense as 
the one by proposed by Ginzburg and Sobyanir£3. Equa- 
tion (56) demonstrates how the critical exponent P — 1/2 
(from (no) oc ~ 1^1^^) can be connected to = 1/3 

(from Ps ~ I^P"^) on account of the comoving density pc. 
The divergent mass coefficient damps the critical fluctu- 
ations such that Eq. (|6|) becomes scaling invariant (this 
is discussed in more detail in Ref. 14). The scaling in- 
variance implies that Eq. (^6|) might be used down to 
|i| = and that the critical exponent of kc and ps might 
be indeed exactly 2/3. This possibility is supported by 
the excellent fits obtained from Eq. (|5J) with Eq. (|5^). 



C. Fit to the experimental data 

We insert the leading t = a' \t\ and kc — ki in Eq. 
i|). This yields 



Pl = a,\t\y^ + a2\t\ + a,\tf^' 
P 



(MAS). (60) 



With 3 parameters we take this as our asymptotic model 
fit, called MAS. Figure^ shows that MAS yields an excek 
lent reproduction of the data by Greywall and AhlersE^ 
(GA), for saturated vapour pressure. 




FIG. 2. Asymptotic model fit of the superfluid density. 
The difi'erences between the 3-parameter fit formula Eq. ( pCj ) 
and the data points are compared to two standard deviations 
(dashed lines). 
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The 3-parameter standard fit (SF) used by GA is 

^ = ai|ir^(l + a3|ir/' + ...) (SF). (61) 

We compare both fits, MAS and SF, by calculating the 
sums of the quadratic deviations: 



XSF 

Xmas 
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(|t| < 0.03) . 



(62) 



This shows that the reproduction shown in Fig. |^ is not 
a trivial result; in the considered range {\t\ < 0.03) the 
data cannot be reproduced by SF. We remark that SF 
(but probably also other 3-parameter fits) fit the data in 
the considerably smaller range \t\ < 0.004; this range is 
used by GA for the fit. The range \t\ < 0.03 seems to 
be appropriate for a 3-parameter fit because already a 
1-parameter fit (ai|iP^^) roughly reproduces the data in 
this range (see Fig. |l|). 

The applicability range of Eq. (jsj) can be considerably 
extended by using the following 4-parameter expansion 
for r and kc'. 

Tit) = a' \tl xS) = ^1 \t\'/' + ^2 |t| + X3 |t|4/3 . 

(63) 

The form of a;c = A kc (t) / ^/In suggests itself because it 
leaves the expansion ( |60| ) unchanged. The 4-parameter 
model fit, Eq. (jsj) with Eq. (|6^), reproduces the data 
down to 1 K within the experimental errors. The fit yields 
the parameter values a' = 3.019, xi — 2.7028, X2 — 
—0.837 and X3 = —3.842; these values have been used 
for Fig. |l|. More details of this fit and a discussion with 
respect to the quasi-particle picture are given in Ref. |lj. 

The fit to the experimental data fixes the decompo- 
sition of ps into po and pc- This decomposition is the 
starting point for the evaluation of the superfluid entropy 



VI. SUPERFLUID ENTROPY 
A. Model prediction 

In Ps = po+ Pc only the condensate part po corresponds 
to a macroscopic wave function and has thus zero entropy 
content. The comoving part pc is made up by different 
s.p.f. (however with the same phase factor in Eq. (^)) 
and has therefore a nonvanishing entropy content. We 
determine this entropy. 

We start with the well-known entropy expression S{nk) 
for a Bose gas with occupation numbers Uk- The IBG 
equilibrium entropy Sibg is obtained from this expres- 
sion by ^iBG = (Sink)) = 5((nfe)iBG). In Sec. |VA 
we have discussed that for r.p.o. in one direction only 
every second s.p. state is occupied. Since the total par- 
ticle number iV = ^ is fixed we have to put twice as 



many atoms in every second state (as compared to the 
plain IBG expression). This leads to the following model 
entropy ^m, 

Sm{T, V,N) = ^Y1 [(1 + 2(^ifc» ln(l + 2{nk)) 



-2{nk) ln(2(nfe)) 

For i = this yields Sm/N = 0.96 fcE 
5ibg/A^ = 1-28 fee. The difference A5(p.o.) : 
has been compared to AE 



(64) 

instead of 
= Sibg — Sm 
< 

!)■ 

in 



[VA. For t 



in Sec 

the (rifc) with T{t) of Eq. ( pSj ) will be used in Eq. (|& 
The resulting overall behaviour of Sm{T) is then — 
contrast to SmciT) — similar to that of the experimen- 
tal entropy S{T). 

The entropy of the superfluid component is Sc — S{q, 
For calculating Sc we have to restrict the sum in Eq, 
by fc < kc, and to account for the c.p.o. (pQ). Prescribing 
the phase field (Pk for each s.p.f. in Eq. ([48| ) implies a two 
to one restriction for each of the three directions of V<?fc . 
For each s.p. state k this reduces the entropy by k^ In 8 
(as long as ^ 1). With these specifications we obtain 

Sc{T, V,N)^'^Y1 [(1 + ln(l + 2(nfe)) 



k<kc 



2{nk) ln(2(nfc))-ln8 . (65) 



The ratio Sc/Sm is the AIBG prediction for the super- 
fiuid entropy fraction Sg/ S, 



Ss Sc 

S Sm 



(66) 



The calculated result is shown in Fig. ||. 



1.8 T(K) 1.7 




FIG. 3. Model prediction for the superfluid entropy Ss as 
a function of the temperature. The figure displays the ratios 
(dashed line) and {S^ / N^) / {S / N) (solid line) both multi- 
plied by a factor 100. The error bars indicate the uncertainty 
of the fit parameters in Eq. (|63|). The prediction becomes 
increasingly uncertain for |t| > 0.1 because the calculation is 
based on an asymptotic expansion 
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B. Experimental detectability 

We show that the AIBG prediction for 5s 7^ is at the 
border of present-day experimental detectabihty. The 
most prominent experiment showing that Ss = (or at 
least Ss « 0) is the fountain pressure (FP) measurement. 
For two containers of liquid helium connected by a su- 
perleak a temperature difference dT produces a pressure 
difference dP. Admitting 5*8 7^ this FP is given by 



dP\ 



V 



1 



S/N 



(67) 



For 5s = this reduces to the well-known London rela- 
tion. The correction term is 5s/5 multiplied by N/Ng = 
p/ Ps] the model prediction for it is shown in Fig. 

The most accurate FP measurements are that by 
Singsaas and AhlersEEl. These authors assume the valid- 
ity of the London relation and interpret their experiment 
as an entropy measurement. For this entropy 5fp, 



V [dTj^p 



(68) 



they find within the errors no deviation from the true 
(caloric) entropy 5. The absolute values of 5 near T\ 
are, however, uncertain by about 2% (Refs. [l^Jl7| ). 

TABLE L Comparison between the caloric (S) and the 
fountain pressure (S'fp) entropy. The input is the experimen- 
tal S'Fp(t) and Cp{t); the quantity AS = S\ — S is calculated 
from Cp. The last column shows the resulting marginal evi- 
dence for S ^ Sfp or, equivalently, for Ss / 0. 



1 


2 


3 


4 


5 


1^1 


pSpp 


Sfp 
Nk^ 


5fp -I- AS 


100 (5 -Sfp) 


Jcm-3K-l 


NkB 


S 


0.0006906 


0.2299 


0.7571 


0.7619 ± 0.0008 


0.07 ±0.11 


0.000790 


0.2296 


0.7561 


0.7615 ± 0.0008 


0.11 ± 0.11 


0.001013 


0.2294 


0.7555 


0.7623 ± 0.0008 


0.02 ± 0.11 


0.001794 


0.2275 


0.7493 


0.7607 ± 0.0009 


0.23 ± 0.11 


0.003177 


0.2250 


0.7411 


0.7601 ± 0.0010 


0.30 ±0.13 


0.005662 


0.2208 


0.7274 


0.7592 ±0.0011 


0.43 ± 0.14 


0.01007 


0.2145 


0.7069 


0.7596 ±0.0013 


0.39 ±0.17 


0.03338 


0.1854 


0.6116 


0.7578 ± 0.0021 


0.74 ± 0.34 


0.07973 


0.1407 


0.4648 


0.7547 ± 0.0034 


1.63 ±0.71 



The theoretical prediction of Fig. || suggests that one 
should compare the temperature dependences of 5 and 
5fp near Tx rather than the absolute values. This is done 
in Table I in a number of steps: 

1. We start with the experimental valueJll for Spp 
in column 1 and 2. By using p{t) from equation 
(Al) of ReL |l8| we relate the entropy to the particle 
number rather than to the density (column 3). 

2. For column 4 we calculated A5 = 5a - S{T) = 
J dT Cp /T from the experimental specific heat Cp 



of Ref. g. If 5fp were equal to 5 then column 
4 should show the temperature independent S\ = 

Sin). 

3. From the temperature dependence of column 4 we 
deduce Sx/N = 0.7624 fee as the limit of 5fp + AS 
for i -> 0". Then 



5-5 



FP 



5a- [5Fp(i)±A5(t)] 5s/7Vs 



5 



AS{t) 



S/N 



(69) 



may be calculated from column 4. The last equality 
in Eq. ^ follows from Eqs. (|^ and (||). Column 
5 displays the resulting experimental evidence (or 
nonevidence) for 5fp 5. 

4. The eriprs included in the Table are the statistical 
errorslij of 0.1% for Swf,{t), and a 1% error for Cp. 
The systematic errorst^ of 5Fp(i) are not included; 
they should be less important because we consider 
only the t-dependence (and not the absolute values) 
in a relatively small interval. 

Figure ^ compares column 5 of the Table with our the- 
oretical prediction. There seems to be some indication of 
a deviation 5 7^ 5fp (t) . The error bars (from the Table) 
show, however, that the experimental evidence for such a 
deviation is at most marginally significant. In any case, 
the compilation of Table I and Fig. |^ shows that, and in 
which way, our theoretical prediction is within the reach 
of experimental detectability. 




FIG. 4. Difference between the caloric (S) and the fountain 
pressure (Sfp) entropy (dots with error bars, from column 5 of 
Table 1). This difference is compared to the model prediction 
(solid line) for {Ss/N,)/(S/N). 



VII. CONCLUDING REMARKS 

We have proposed a microscopic model for the A-tran- 
sition of liquid helium. The model uses phenomenologi- 
cal assumptions, in particular the transition itself is in- 
troduced by the analogy to the IBG. The basic idea of 
the AIBG is that phase ordering leads to an extra en- 
ergy ^ J2 which — evaluated with the IBG occupa- 
tion numbers — yields a logarithmic singularity. In this 
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model the superfluid density is not identical to the square 
of the order parameter field {ps ^ rn\'4'\'^)- A fit to the 
experimental ps leads to the model prediction for a non- 
vanishing superfluid entropy Ss- 

The model is hardly related to standard models of liq- 
uid helium which are either microscopic variational ap- 
proaches (like Rcfs. ||,^^Lor the Landau- Wilson renor- 
malization group theoryE3. Because of its novelty and 
originality our model cannot compete in quality of foun- 
dation, exactness and completeness with those other the- 
ories. It leads, however, to specific predictions about the 
structure of ps. These predictions are without compe- 
tition by other approaches and they can and should be 
tested experimentally. Therefore the presentation of the 
model seems appropriate in spite of a number of unsolved 
questions concerning the model assumptions. 

In contrast to other approaches the AIBG offers also 
a possible solution of the so-called microscopic prob- 
lem of liquid helium. This problem was formulated by 
UhlenbeckEil: If ps is identified with a single quantum 
state (ps = to|?/'P) then the approach to equilibrium 
(Ps = Ps{T)) cannot be understood. In the AIBG this 
problem is solved by the contribution of noncondensed 
particles to the superfluid density. 
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